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ABSTRACT 


The  analysis  of  a  plane  wave  propagating  through  a  homogeneous 
plasma  region  of  dimensions  large  compared  to  the  signal  wavelength 
has  shown  that  the  wave  attenuation  can  be  related  to  physical  param¬ 
eters  of  the  plasma  such  as  electron  density,  temperature,  and  pres¬ 
sure.  In  this  report,  the  plasma  in  the  combustion  chamber  of  a  small 
solid  propellant  rockec  motor  is  used  to  investigate  the  feasibility 
of  determining  the  temperature  of  the  plasma  by  using  microwave  atten¬ 
uation  data. 

Potassium  perchlorate  is  substituted  for  some  of  the  ammonium 
perchlorate  in  the  propellant  to  give  an  electron  decsity  that  pro¬ 
duces  an  easily  measurable  attenuation.  The  theoretical  results  show 
that  for  a  known  propellant  composition  and  potassium  perchlorate 
concentration,  both  the  plasma  temperature  and  chamber  pressure  at 
equilibrium  may  be  found  from  one  attenuation  measurement.  Also,  for 
a  known  t  ;ermodynamic  state  of  the  plasma,  the  electron  density  may 
be  found  ft  fa  the  attenuation  data. 

The  method  of  microwave  analysis  is  shown  to  be  very  useful  in 
the  study  of  transient  phenomena  in  that  nc  transducers  are  required; 
hence,  the  time  response  of  the  measurement  is  limited  only  by  the 
microwave  detection  circuitry.  This  technique  may  be  used  to  study 
localized  disturbances  within  the  plasma  region. 

Experimental  data  are  presented  to  support  the  theory  for  an 
ammonium  perchlorate-polybutadiene  acrylic  acid  propellant  which  was 


burned  at  chamber  pressures  between  20  and  80  psia  corresponding  to 
a  temperature  range  of  2550°K  to  2600#K.  The  data  show  that  temper¬ 
ature  changes  of  les£<  than  ten  degrees  can  easily  be  detected  by 
microwave  attenuation  measurements  wher.  the  propellant  is  seeded  with 
2-10  percent  potassium  perchlorate.  Seedjng  other  propellants  with 
alkali  metal  compounds  is  expected  to  yield  comparable  results  over 
the  combustion  range  in  temperature  and  pressure  of  the  other  propel¬ 


lants. 
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I.  DESCRIPTION  OF  THE  INVESTIGATION 


1.1  Introduction 

The  difficulty  of  measuring  the  temperature  of  a  hot  plasma 
by  conventional  probe  transducers  has  indicated  the  need  for  another 
temperature  measuring  method.  A  method  that  uses  microwave  tech¬ 
niques  to  analyze  plasma  parameters  has  been  developed1  This  method 
has  been  applied  to  the  analysis  of  the  exhaust  plume  of  a  rocket  and 
yields  information  about  the  plasma  temperature,  pressure,  and  elec¬ 
tron  density  from  changes  in  the  propagation  constants  of  an  incident 
electromagnetic  wave.  The  present  investigation  concerns  a  modifica¬ 
tion  of  this  method  and  considers  a  small,  solid  propellant  rocket 
motor  in  which  quartz  windows  are  used  to  allow  passage  through  the 
combustion  chamber  of  a  focused  microwave  signal. 

The  propellant  is  seeded  with  various  amounts  of  potassium 
perchlorate  to  furnish  a  significant  electron  density  to  enhance  the 
electromagnetic  wave  interaction.  It  is  assumed  that  the  partially 
ionized  potassium  is  the  only  source  of  free  electrons,  and  the  temper¬ 
ature  of  these  electrons  can  be  determined  from  a  measurement  of  the 
wave  attenuation  and  a  knowledge  of  the  chemical  composition  of  the 
propellant. 

The  theoretical  basis  for  the  determination  of  the  plasma  param¬ 
eters;  temperature,  pressure,  and  electron  density  is  given  in  the 

1  A.  S.  Jones,  C.  C.  Johnson,  and  R.  W.  Grow,  IEEE  Transactions,  Nuclear 
Science,  Vol.  NS-11,  January  1964,  p.  277. 


following  section.  The  experimental  apparatus  is  described  in  Chapter 
IZ  and  the  results  and  conclusions  in  Chapter  III.  Zt  will  be  seen 
that  the  presented  analysis  gives  instantaneous  information  about  the 
plasna  and  is  veil  adapted  to  the  study  of  conditions  in  various  re¬ 
gions  of  the  chamber. 

1.2  Theoretical  Basis 

The  interaction  of  an  electromagnetic  wave  with  a  high  tempera¬ 
ture,  multispecies  plasma  is  discussed  in  this  section  In  terms  of  the 
attenuation  of  an  incident  plane  wave  propagating  through  the  combus¬ 
tion  chamber  of  a  solid  propellant  rocket  motor. 

The  dispersion  equation  describing  the  propagation  of  the  wave 
through  the  plasma  region  is  first  derived  using  Maxwell's  equations 
and  the  Langevin  force  equation.  These  mathematical  relations  will 
describe  the  wave  attenuation  in  terms  of  the  frequency  of  the  wave, 
the  electron  density  of  tbe  plasma,  and  the  electron  drift  velocity. 
The  physical  parameters  of  the  drift  velocity  will  be  shown  to  give 
the  mean  electron  velocity  as  a  function  of  the  species  of  gas  present 
and  their  electron  collision  cross  sections,  gas  pressure  and  tempera¬ 
ture,  and  wave  frequency. 

In  the  propellant,  potassium  perchlorate  was  substituted  for 
some  of  the  ammonium  perchlorate.  Pot  ass  list,  having  a  relatively  lo» 
ionization  potential,  furnishes  enough  electrons  to  swamp  the  random 
generation  of  free  electrons  and  produce  a  measurable  attenuation  in 
the  pressure  range  of  20  ?sia  to  100  psia.  The  final  form  of  the  at- 
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That  Eq.  1.4  is  a  scalar  relation  can  be  assumed  from  the  geometry 
of  the  test  motor  and  the  frequency  of  the  applied  signal.  The 
quarts  window  apertures  of  the  combustion  chamber  measured  one  by 
two  inches;  and  since  the  22.275  GHz  frequency  wave,  having  a  free 
space  wavelength  of  0.519  inches,  was  focused  using  wave  guide  horns 
and  dielectric  lenses,  the  plane  wave  analysis  appears  to  be  well 
justified. 

A  damping  term  is  used  in  the  Langevin  force  equation  to 
account  for  the  energy  loss  of  the  electron-particle  collisions.3  The 
force  equation  is 

v  +  gv  «  e/a  E  (1.5) 

which,  for  a  given  harmonic  time  dependence,  has  a  solution  of 

V  *  e/*  g'fg  (1*6) 

For  an  electric  field  intensity  having  a  simple  one-dimensional 
space  dependence  given  by  . 

E  »  E  e"rz  (1.7) 

o 

the  simultaneous  solution  of  the  Langevin  equation,  Eq.  1.6,  the  wave 

3  Ibidt  p.  394. 
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!  equation,  F-q.  1.4,  and  Eq.  1.7  for  the  propagation  constant,  r,  yields 

|  the  dispersion  equation. 


r2  +  k2  -  k2 


pg  +  i» 


(1.8) 


where  k  is  the  propagation  constant  of  free  space  and  k^  is  defined 
by 
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k  -  -E 
P  c 


(1.9) 
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The  quantity  u 


is  the  plasma  frequency  squared  and  is  given  by 


2 

P 


(1.10) 


the  quantity,  p£,  is  the  electron  charge  density  of  the  medium  in 
coulombs  per  unit  volume.  In  the  dispersion  equation,  the  complex 
quantity  f  is  usually  separated  such  that 


r  -  u  +  jB 


(1.11) 


where  a  is  the  attenuation  constant  in  nepers  per  meter  and  g  is  the 
phase  constant  in  radians  per  meter. 

The  parameter  g  in  the  *orce  equation  has  dimensions  of  fre¬ 
quency,  and  for  a  simple  plasma  where  g  is  independent  of  electron 
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velocity,  it  reduces  to  the  collision  frequency  of  the  electrons.4  The 
combustion  chamber  plasma  assumed  constitutes  a  very  hot,  dense  mixture 
of  gases,  and  under  these  conditions,  the  quantity  g  may  be  complex  and 
unequal  to  the  collision  frequency.  In  order  to  allow  for  this  possi¬ 
bility,  the  solution  to  the  force  equation  may  be  written  as 

v  -  —  (B  -  iD)  (1.12) 

m 

a  fora  used  by  Mo lmud. 5  The  quantities  B  and  D  can  be  related  to  the 
conductivity  of  the  plasma  by  noting  the  expression  defining  current 
density  and  combining  this  relation  with  Eq.  1.3,  Eq.  1.10,  and  Eq. 

1.12.  The  current  density  is 


J  »  oE 


(1.13) 


And  the  use  of  this  equation  with  the  other  relations  yields  the  medium 
conductivity. 


a  “  (B  -  iD) 

p  o 


(1.14) 


4  M.  A.  liman.  Introduction  to  llama  Physics ,  McGraw-Hill  Book  Company, 

New  York,  1964,  p.  42. 

5  p.  Mo  lmud ,  "Langevin  Equation  and  A  C  Conductivity  of  Hon-Maxwellian 

Plasmas, "  Physical  Review,  Vol.  114,  April  1959,  pp.  29-32. 
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The  modified  dispersion  equation  is  then 


2  2  2  2 
r  +  1c  -  jk  mB  -  k  uD 
P  P 


(1.15) 


If  the  propagation  constant  T  is  represented  as  in  Eq.  1.11,  the  de~ 
terminental  equation  of  Eq.  1.15  can  be  separated  into  real  and  imag¬ 
inary  parts: 
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1/2 


1/2 


(1.16) 


(1.17) 
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It  is  seen  that  a  and  3  are  functions  of  the  frequency  of  the  wave, 
the  plasma  frequency  which  is  related  to  electron  density,  and  the 
complex  drift  velocity  terms,  B  and  D. 

1.4  Investigation  of  the  Drift  Velocity 

In  order  to  present  the  attenuation  and  phase  constants  in 
terms  of  the  measurable  parameters  of  the  system,  the  functional  de- 
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where  fQ,  f ^ ,  and  are  coefficients  of  the  expansion.7  Equation 
1.18  may  then  be  written  as 


,2it  *  • 


III  vz[(fi  cos  Ut  +  Ei  sin  wt)  vzj  v^dv  sin6  d8  d+ 


4  -  0  6*0  v  “0 
z 


2it  *  • 

I  I  I 

4*0  0=0  v  *0 
z 


(1.21) 


f  v  dv  sin0  d8  d4 
o  T 


In  addition,  relations  may  be  developed  among  f  ,  f,,  and  g,  such  that 

o  1  1 

Eq.  1.21  may  be  expressed  in  complex  form  as 


3/2  _  mv_ 
2kT 


1 1 1  *' '**  7777' fe)  •'  ”  (>-») 


cos  0  sine  d0  d+ 


Hi-fi 


.3/2  av 

2St)  e  2kT  v2  8ine  d0  dv  d* 


(1.22) 


where  X  is  the  mean  free  path  between  collisions  and  n  is  the  electron 
density.8  If  a  change  of  variable  is  made  such  that 


7  H.  Margenau,  "Conduction  and  Dispersion  of  Ionized  Gases  at  High  Fre¬ 

quency,”  Physical  Review,  Vol.  69,  May  1946,  pp.  575-585. 

8  Ibid. 


-  9  - 


(1.23) 


v 


v 

X 


and  the  8  and  $  integrations  are  carried  out,  Eq.  1.22  becomes 


where 


v 

z 


ff /  • 


f/2  u4  e'^’2  dv 
J  i&i  +  V 


s 


m 

2kT 


(1.24) 


(1.25) 


For  so-called  Maxwellian  gases,  in  which  the  potential  between 

-4 

electrons  and  neutral  molecules  varies  as  r  ,  the  collision  fre¬ 
quency  is  independent  of  7,  allowing  inmediate  integrations  of  Eq.  1.24; 
however,  for  a  mixture  of  various  types  of  gases,  the  average  colli¬ 
sion  frequency  is  found  from  the  average  electron-molecule  collision 
cross  section  given  by 


Q  ■  7  m.Q.  (1.26) 
avg  L  ixi 

where  the  's  ate  the  mole  fraction  averaging  factors  of  each  species 
of  gas  present.  Since  the  collision  frequency  of  the  electrons  and  a 
particular  neutral  gas  is  given  by 
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vi  “  pvQi 


(1.27) 


where  p  is  the  gas  particle  density.  The  average  collision  frequency 
for  a  given  electron  velocity  is 


v  »  pvQ 
avg  K  avg 


(1.28) 


Substitution  of  Eq.  1.26  in  Eq.  1.28  gives 


v  *  pv  V  m.Q. 
avg  K  *•  ixi 


(1.29) 


The  average  drift  velocity  of  electrons  due  to  the  incident 
microwave  signal  as  given  in  Eq.  1.24  is  now  a  function  of  the  fre¬ 
quency  of  the  wave,  the  mole  fraction  of  each  gas  species,  and  each 
collision  cross  section: 


H  -p*  , 

v  e  dv 
ia>  +  pv  l  »1Q1 


(1.30) 


To  obtain  useable  expressions  for  the  cross  section  terms  in 
the  summation  in  Eq.  1.30,  a  power  series  approximation  of  the  colli¬ 
sion  cross  sections  will  be  used.  This  approximation  is  valid  for 
only  a  limited  range  of  velocity;  however,  it  is  known  that  the  varia- 
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tion  in  the  combustion  gas  temperature  for  pressures  between  20  psia 
and  100  psia  is  only  about  50°K  out  of  an  equilibrium  temperature  of 
2600*K,  indicating  that  the  thermal  energy  of  the  electrons  actually 
changes  very  little  over  the  entire  pressure  range  considered.  It 
will  be  shown  in  Chapter  II  that  the  significant  collision  species 
in  the  plasma  are  water,  hydrogen  chloride,  carbon  dioxide,  carbon 
monoxide,  diatomic  hydrogen,  and  diatomic  nitrogen.  Abundant  informa¬ 
tion  exists  about  the  dependence  of  the  collision  cross  section  on 
electron  energy  for  these  substances.9*10*11  Using  the  fact  that  the 
thermal  energy  of  the  plasma  electrons  is  of  the  order  of  0.22  elec¬ 
tron  volts  and  varies  only  a  few  percent  from  that  value,  the  follow¬ 
ing  expressions  represent  good  approximations  to  the  collision  cross 
sections  for  the  given  gaseous  species: 


5.30  x  1Q~8 
H2°  *  v2 


(1.31) 


n  »  1*835  x  IQ 

%C1  "  2 

v 


(1.32) 


9  I.  P.  Shkarofsky,  T.  W.  Johnson,  and  M.  P.  Bachynski,  The  Particle 

Kinetics  of  Plasmas ,  Addison-Wesley Publishing  Company,  Inc., 
Beading,  Pennsylvania,  1966,  pp.  175-184. 

10  H.  S.  W.  Massey  and  E.  H.  S.  Burhop,  Electronic  and  Ionic  Impact 

Phenomena ,  Oxford  University  Press,  London,  England,  1952,  pp. 
205-210. 

11  S.  Altschuler,  "Theory  of  Low-Energy  Electron  Scattering  by  Polar 

Molecules,"  Physical  Review,  Vol.  107,  No.  1,  July  1957,  pp. 
114-117. 


-  12  - 


(1.33) 


-1/2 

y  -  8  v 


(1.39) 


and  substituting  Eq.  1.39  and  Eq.  1.12  into  Eq.  1.38  gives 


B-iD 


8  f  8y3  e  y  dy  _ _ 

3*^  i«8y  +q2c2p®^^  +  mic^P8y  +mocop8ii,2y2  +  m_1c_1py3 


(1.40) 


Upon  rationalizing  the  denominator ,  it  is  seen  that 


•  (  c3/2  a.  o  *1/2  2  .  3V,  5  -y  . 

8  f  Vm2C2pg  +Ylp8ytB0Vf  y  +n-ic-lpy  jfiy  e  dy 

1*^*  J  2  /  3/2  1/2  2  3\^ 

o  («8y)  +^»2c2pB  +a1c1p8y+«ocopB‘/  y  +m_1c_1py  | 


(1.41) 


|  _ »82y^  e  ^  dy _ 

o  («8y)2+r«2c2p83^2+m1c1p8y  +  nocop81^2y2  +  m^c^py3^ 


(1.42) 


Uith  the  determination  of  B  and  D  as  functions  of  the  incident 
wave  frequency,  gas  composition,  and  temperature,  the  propagation  con¬ 
stants  of  Eq.  1.16  and  Eq.  1.17  are  determined  except  for  the  value  of 
the  plasma  frequency.  Once  the  electron  density  is  known  for  the  plasaa 
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of  a  given  propellant,  the  'Jave  attenuation  is  completely  determined 
and  can  be  theoretically  predicted  for  a  given  geometry. 


1.5  Investigation  of  the  Electron  Density 


In  this  derivation  of  the  propagation  constants  for  the  inci¬ 
dent  electromagnetic  wave,  the  temperature  contained  in  the  expres¬ 
sions  thus  far  is  the  temperature  of  the  electrons  and  was  introduced 
through  the  moan  velocity  expression,  Eq.  1.12.  When  the  electrons 
are  in  thermal  equilibrium  with  the  gas  molecules,  the  temperature  in 
a  and  8  can  be  equated  to  the  gas  temperature,  which,  in  turn  relates 
tc  the  gas  density  through  the  perfect  gas  law 


(1.43) 


where  Nq  is  Avogodro’s  number,  R  is  the  gas  constant,  and  p  is  the 
pressure.  The  perfect  gas  law  then  allows  the  expression  of  density 
in  terms  of  temperature  provided  the  pressure  is  known.  The  equilib¬ 
rium  electron  density  then  may  be  derived  from  the  thermal  state  of 
the  gas.  It  will  be  assumed  that  the  source  of  electrons  is  the  seeded 
potassium  perchlorate  and  that  the  mechanism  of  free  electron  produc- 
ticn  is  thermal  ionization.  There  are  two  ways  of  considering  the  pro¬ 
duction  of  electrons;  they  are: 

1.  Assuming  that  no  sinks  for  electrons  exist  in  the  gas. 

2.  Assuming  that  chlorine  acts  as  a  sink  for  electrons. 
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Assuming  No  Electron  Sinks 


For  the  simple  reaction 


X  -►  K+  +  e" 


(1.44) 


Saha's  equation  may  be  written  as 


l  -  x 


(1.45) 


where  N  is  the  total  number  density  per  unit  volume  of  potassium* 

X  is  the  fraction  of  those  atoms  ionized*  k  is  Boltzmann's  constant* 
h  is  Plank's  constr^t,  is  the  electronic  mass*  e  is  the  electronic 
charge*  and  I  is  the  ionization  potential.12  For  a  given  temperature, 
the  number  density  for  electrons  per  unit  volume  is  simply 


n  -  XH 
e 


(1.46) 


and  the  plasma  frequency  defined  in  Eq.  1.10  becomes 


2 

u 

P 


a  e 
e  o 


(1.47) 


12  S.  Altschuler*  P.  Mo laud ,  and  M.  Hoe,  "The  Electromagnetics  of  the 
Rocket  Exhaust*'1  Space  Technology  Laboratory  Report,  GM  TR  0165- 
00397,  June  1958,  p.  130. 


Consideration  of  Electron  Sinks 


When  electron  recombination  is  taken  into  account,  not  only  the 
number  density  of  the  electron  source,  but  also  the  number  density  of 
the  electron  sink  must  be  known.  Chlorine  will  be  assumed  to  be  the 
only  effective  electron  sink  because  it  is  the  only  constituent  pres¬ 
ent  in  substantial  amounts  in  the  propellant  that  has  an  electron  af¬ 
finity  of  the  order  of  the  low  ionization  potential  of  potassium. 

With  the  following  two  reactions  considered,  the  chlorine  elec¬ 
tron  sink  is  defined: 


K  -*•  K+  +  e 


(1.48a) 


Cl  +  e  -►  Cl 


(1.48b) 


The  chemical  equilibrium  constants  for  this  pair  are  given  by 


[did 

IKJ 


(1.49a) 


[ci] 

U-] 

[< 

Cl'] 

(1.49b) 


When  the  bracketed  quantities  denote  concentrations  in  number  per 
unit  volume,  the  equilibrium  constants  may  be  written  as 


and 


(l.50a) 


(1.50b) 


where  1^  is  the  ionization  potential  of  potassiua  and  is  the 
electron  affinity  of  chlorine.13 

Defining  the  total  quantities  of  the  reacting  elements  by  the 
following  equations: 


A1-[k]+[r+]  (1.51a) 

and 

*2  «  [cl]  +  [ci“]  (1.51b) 

it  is  possible  to  coabine  Eqs.  1.49  and  1.51  with  the  neutrality  con¬ 
dition. 


H +  N  -  [*+j 


(1.52) 


13  S.  Glssstone,  Textbook  of  Physical  Chemistry,  Van  Nostrand  Coapany, 
New  York,  19*7,  p.  882. 


-  18  - 


Co  tom  a  cubic  equation  in  e.ectron  density, 

n^  +  An^  +  Bn  +  C  *  0 


(1.53) 


where 


*  *  M 


(1.54) 


and  the  coefficients  are  given  by 

A  *  Oj  +  o2  +  Aj  (1.55a) 

•  ■  -iK  -  *i +  “2) 

C  “  -OjC^Aj 

The  determined  value  then  establishes  the  plasma  frequency  through 
Eq.  1.4C. 

1 .6  Plasma  Temperature  frog 
the  Attenuation  Constant 

To  relate  the  attenuation  of  a  microwave  beam  in  a  plasma  to 
the  temperature,  consider  that  a  plane  wave  normally  impinges  upon 
a  plasma  region  and  passes  through  the  plasma  into  free  space  as 
illustrated  in  Fig.  1.  At  plane  I  in  the  figure,  due  to  differences 
in  the  characteristic  impedances  of  the  plasma  and  free  space,  a  re- 


(1.55b) 

(1.55c) 
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Fig.  X.  Diagram  of  tha  aaaumad  alactromagnaeic  fiald  gaomaery. 


flection  coefficient  nay  be  defined  in  terns  of  the  two  impedances  by 


R 


Z.  -  Z 
1  o 

Z.  +  Z 
1  o 


(1.56) 


where  Zq  is  the  characteristic  impedance  of  free  space  and  is  that 
of  the  plasma.  If  it  is  assumed  that  the  reflection  coefficient  is 
small  compared  to  unity,  then  the  amplitude  of  the  transmitted  wave 
an  infinitesimal  distance  to  the  right  of  plane  I  will  be  essentially 
equal  to  the  amplitude  of  the  free  space  wave.  If  the  transmitted 
amplitude  of  the  electric  field  is  designated  as  and  the  incident 
wave  by  E^,  then  the  amplitude  of  the  electric  field  within  the  plasma 
is  given  by 


E  -  E  e~r*  (1.57) 

p  o 

where  z  is  measured  to  the  right  of  the  plasma  boundary  I.  If  the 
assumption  of  a  small  reflection  coefficient  is  again  made  at  plane 
II,  then  the  electric  field  amplitude  at  plane  II  designated  by  Ej 
is  given  by 


*r 


E  e 
o 


-n 


(1.58) 


The  attenuation  of  the  transmitted  wave  in  decibels  by  definition  is 


-  21  - 


£ 

£ 

jr 

f5 


or 


20  log1()  f- 


(1.59a) 


A  -  20  log10  c 


-fL 


(1.59b) 


In  tens  of  the  magnitude  of  a  of  Eq,  1.16, 


A  *  20  aL  log^Q  e 


(1.60) 


Substitution  of  Eq.  1.16  into  Eq.  1.60  yields 


A 


1C/T  kL  log10 


It  is  seen  that  one  attenuation  measurement  through  a  specified  path 
length  of  plasma  of  a  known  cosposition  and  pressure  can  yield  the 
plassa  temperature. 
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II.  APPLICATION  OF  THE  THEORY 


In  chis  chapter,  the  relations  derived  in  Chapter  I  are  applied  to  a 
specific  test  situation.  Inasmuch  as  the  theoretical  attenuation  de¬ 
pends  on  the  geometry  of  the  combustion  chaaber  and  composition  of  the 
propellant,  the  test  unit  will  be  described  first.  Second,  a  descrip¬ 
tion  will  be  given  of  the  computer  program  that  was  used  to  solve  for 
the  theoretical  wave  attenuation.  Finally,  the  results  of  the  analy¬ 
sis  will  be  presented. 

2.1  The  Test  Unit 

Figure  2  shows  two  photographs  of  the  test  rocket  motor  used  in 
this  investigation.  The  inside  diameter  of  the  combustion  chamber  was 
two  inches  while  the  perpendicular  distance  between  the  two  quartz 
windows  measures  2.20  inches.  Various  nozzle  diameters  were  used  to 
produce  pressures  ranging  from  approximately  80  psig  for  a  0.235  inch 
diameter  nozzle  to  18  psig  for  a  0.3125  inch  diameter. 

The  solid  propellant  used  in  the  test  motor  was  a  variation  of 
ammonium  perchlorate  (AP)  and  polybutadiene  acrylic  acid  (PBAA)  bonded 
with  Epon  826  resin,  the  percentage  by  weight  given  by  82  percent,  15.3 
percent,  and  2.7  percent,  respectively,  and  the  variation  being  a  sub¬ 
stitution  of  potassium  perchlorate  ^KCIO^)  for  some  of  the  AP.  Four 
concentrations  of  potassium  perchlorate  were  used;  1  percent,  2  percent, 
5  percent,  and  10  percent  weigtst. 

Pressure  war  measured  with  two  devices;  a  Kistler  gauge  was  in¬ 
serted  into  the.  rear  of  the  chamber,  and  a  Statham  gauge  was  connected 
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to  the  port  above  the  windows  as  shown  in  the  photographs. 

The  microwave  portion  of  the  test  unit  consisted  of  a  klystron 
oscillator  wave  guide  fed  to  the  chamber  through  a  1  kHz  ferrite  modu¬ 
lator.  The  wave  guide  horns  and  dielectric  lenses  used  for  focusing 
the  beam  can  be  seen  in  the  photograph  (a).  After  passage  through  the 
coabusion  chamber,  the  K-band  signal  was  sent  through  a  series  of  at¬ 
tenuators  to  lower  the  power  level  for  crystal  detection.  One  attenua¬ 
tor  was  a  precision  unit  that  served  as  a  calibration  reference.  The 
detected  signal  was  connected  to  a  power  indicator  and  recorded  on  a 
Sanborn  strip  recorder  along  with  the  Statham  gauge  pressure  reading. 

2.2  The  Computer  Program 

The  solution  of  the  plasma  attenuation  constant  was  carried  out 
on  the  Univac  110b  computer  at  the  University  of  Utah  Computer  Center. 
The  program  listing  is  given  in  the  Appendix,  and  it  is  divided  into 
three  major  parts: 

1.  Initialization  of  constants  representing  a  specific  condi¬ 
tion  of  gas  and  electron  temperature. 

2.  Calculation  of  the  mean  velocity  terms  B  and  D,  electron 
density,  and  attenuation. 

3.  Printing  the  calculated  attenuation  and  other  pertinent 
parameters  of  the  plasma. 

Initialization  of  Constants 

The  initialization  section  defines  physical  constants  such  as 
frequency,  velocity  of  light,  etc.,  and  then  reads  from  the  data  deck 
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the  mole  fraction  concentrations  of  the  important  constituents  of  the 
combustion.  Values  used  for  the  representative  mole  fractions  were 
taken  from  a  free  energy  analysis  of  the  combustion.  The  particular 
method  used  was  a  theoretical  specific  impulse  (ISP)  computer  program 
furnished  to  the  University  of  Utah  by  the  United  States  Air  Force 
Rocket  Propulsion  Laboratory  at  Edwards  Air  Force  Base,  California. 
Given  the  chemical  composition  of  the  propellant  and  a  stated  combus¬ 
tion  pressure,  this  program  will  print  all  the  equilibrium  constants 
and  chemical  components  of  the  resulting  gas.  For  the  propellant  com¬ 
positions  and  chamber  pressures  described  earlier  in  this  chapter,  more 
than  ninety-eight  percent  of  the  equilibrium  mass  is  accounted  for  by 
six  gaseous  species: 

1.  diatomic  hydrogen 

2.  water 

3.  diatomic  nitrogen 

4.  carbon  monoxide 

5.  carbon  dioxide 

6.  hydrogen  chloride 

The  mole  fractions  for  these  constituents  were  taken  from  the  program 
print  out  end  constitute  the  values  in  the  summation  of  Eq.  1.29,  which 
expression  is  used  in  the  calculation  of  the  mean  velocity  terms  S  and 
1).  A  table  of  mole  fraction  values  that  constituted  the  data  deck  is 
also  given  in  the  Appendix. 

After  the  data  are  read  in,  the  temperature  is  initialized  at 
2000  degrees  absolute.  This  temperature  does  not  represent  the  equili- 
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brium  temperature  of  the  gas  predicted  by  the  free  energy  analysis,  but 
it  merely  serves  as  a  starting  reference  temperature  for  attenuation 
calculations.  Duriug  steady-state  conditions  of  combustion,  lw  can  be 
assumed  that  the  eleccrons  and  gas  are  in  thermal  equilibrium,  yet  it 
is  conceivable  that  in  a  transient  situation,  the  temperature  could 
change  from  the  equilibrium  value  and  alter  the  attenuation.  For  this 
reason,  a  range  of  temperatures  in  a  neighborhood  about  the  equilibrium 
gas  temperature  is  used;  hence,  for  a  given  gas  state,  a  number  of  at¬ 
tenuations  will  be  calculated  for  a  group  of  surrounding  temperatures. 
Two  other  chemical  concentrations  are  read  from  the  data  deck.  These 
are  the  equilibrium  amounts  of  potassium  and  chlorine,  the  assumed 
source  and  sink  of  electrons,  respectively. 

Calculation  of  Attenuation  Quantities 

After  the  physical  constants  of  the  plasma  have  been  defined,  it 
is  necessary  to  calculate  three  quantities: 

1 .  The  mean  velocity  terms  B  and  D. 

2.  The  number  density  of  electrons. 

3. _  The  theoretical  attenuation  itself. 

First  is  the  evaluation  of  the  integrals  of  B  and  D  from  Eqs.  1.41  and 
1.42.  The  integrands  are  continuous,  positive  functions,  and  the  in¬ 
tegrals  are  easily  evaluated  by  a  simple  trapezoidal  method. 

The  number  density  of  free  electrons  ia  next  determined  from  Eq. 
1.53  using  a  subroutine  to  solve  for  the  roots  of  the  cubic  equation. 

It  must  be  stated  that  for  some  outlying  temperatures,  the  coefficients 


of  the  coble  equation  were  such  Lhst  round  off  error  in  the  cubic  cal¬ 
culations  gave  no  positive  real  roots.  This  fact  is  equivalent  to  the 
obviously  nonphysical  situation  of  assuming  a  temperature  far  away  from 
the  equilibrium  value.  The  relatively  large  amount  of  chlorine  com¬ 
pared  to  potassium  rendered  the  cubic  equation  root  locus  sensitive  to 
lower  than  equilibrium  tenmeratures .  The  occurrence  of  no  positive 
solutions  was  not  felt  Important  because  the  temperatures  giving  no 
roots  were  always  more  than  one  hundred  degrees  from  the  indicated 
equilibrium  value. 

Earl;  in  the  experiments  of  this  investigation,  it  was  seen  that 
the  assumption  of  no  electron  sinks  gave  predicted  attenuations  that 
were  much  larger  than  those  observed  in  the  test  motor;  therefore,  the 
use  of  Saha's  equation  will  be  assumed  incorrect,  and  the  computer  pro¬ 
gram  uses  only  the  cubic  expression  to  determine  the  electron  density. 

The  Computer  Print  Out 

For  an  assumed  gas  composition  and  temperature  range,  nine  plasma 
quantities  are  printed.  They  are  temperature  in  degrees  absolute,  gas 
density  in  number  per  cubic  meter,  charge  density  in  coulumbs  per  cubic 
meter,  B  and  D  integral  evaluations  from  Eqs.  1.41  and  1.42  in  units 
of  inverse  frequency  or  seconds,  electron  number  density  in  number  per 
cubic  meter,  the  plasma  frequency  as  defined  in  Eq.  1.10  in  units  of 
frequency  or  inverse  seconds,  and  the  attenuation  listed  in  both  nepers 
per  meter  and  decibels  loss  p*;r  a  2.20  inch  attenuation  path. 

Each  printed  page  contains  a  temperature  range  of  2000  degrees 


to  3000  degrees,  and  a  total  of  twenty-four  pages  Is  printed  accouat- 
ing  for  the  four  potassium  perchlorate  concentrations  and  six  chamber 
pressures,  120-,  100-,  SO-,  60-,  40-,  and  20-psia  for  each  KCIO^  seed¬ 
ing.  From  the  computer  print  out,  the  theoretical  results  follow 
directly. 

2.3  Theoretical  Results 

From  the  calculated  attenuations  of  the  computer  program  and 
the  equilibrium  temperatures  predicted  by  the  Air  Force  ISP  program, 
three  plasma  attenuation  characteristics  may  be  derived.  One  set  of 
curves  represents  non-equilibrium  states  of  the  plasma  and  is  a  graph 
of  attenuation  versus  temperature  over  a  temperature  range  about  equi¬ 
librium.  The  other  two  graphs  are  for  an  assumed  equilibrium  state; 
one  shows  attenuation  versus  temperature  for  the  various  compositions 
of  propellant,  and  the  other  shows  attenuation  versus  pressure  for 
different  compositions. 

Figure  3  is  a  plot  of  the  non-equilibrium  states  of  the  plasma 
taken  directly  from  the  values  printed  by  the  attenuation  program  for 
a  KCIO^  concentration  of  10  percent.  Figure  4  is  a  similar  graph  for 
a  2  percent  seeding  of  KCIO^.  It  is  seen  that  the  predicted  attenua¬ 
tion  is  a  strong  functioa  of  both  the  initial  KCIO^  content  and  pres¬ 
sure. 

Figure  b  shows  the  theoretical  attenuation  versus  temperature 
for  the  equilibrium  condition.  From  the  lines  of  constant  pressure 
that  are  shown,  it  is  again  seen  that  the  attenuation  is  a  strong  func- 
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tion  of  pressure.  Figure  6  gives  the  attenuation  versus  equilibrium 
chamber  pressure.  This  graph  is  merely  a  modification  of  Fig.  5  to 
show  more  clearly  the  attenuation-pressure  relationships  since  these 
quantities  were  used  as  measured  parameters  in  the  experimental  tests. 

The  presented  graphs  clearlv  show  how  much  attenuation  is  to 
be  expected  in  a  K-band  signal  incident  upon  a  2.20  Inch  cross  sec¬ 
tion  of  homogeneous  plasma  of  assumed  composition.  The  experimental 
verification  of  these  predictions  and  the  conclusions  of  this  investi¬ 
gation  are  given  in  the  following  chapter. 


III.  EXPERIMENTAL  RESULTS  AND  CONCLUSIONS 


In  this  chapter,  experimental  data  are  given  to  support  the 
theoretical  results  given  in  Chapter  II.  It  will  be  shown  that  for 
the  equilibrium  conditions  of  the  combustion  chamber  assumed  in  the 
theoretical  ISP  program,  the  measured  wave  attenuation  of  the  test 
rocket  motor  closely  agrees  with  the  predictions  of  Figs.  5  and  6. 

Also  given  in  this  section  is  a  brief  summary  of  the  practical  con¬ 
sequences  of  the  results  of  the  investigation  along  with  some  general 
conclusions  of  the  technique  of  using  microwave  methods  in  the  analy¬ 
sis  of  a  combustion  chamber  plasma. 

3.1  Experimental  Verification 

The  test  configuration  has  been  described  in  Chapter  II,  and 
only  the  experimental  results  are  presented  here.  Figure  7  shows  a 
representation  of  the  strip  chart  recording  of  a  typical  firing  of 
the  test  motor.  The  basic  features  are  obvious  from  an  examination 
of  the  figure.  The  slow  shift  in  pressure  has  been  attributed  in  some 
measure  to  acoustic  waves  in  the  chamber;  however,  the  experiment  was 
not  designed  to  analyze  such  high  frequency  pressure  phenomena,  so  the 
slow  changes  in  pressure  during  the  burn  do  not  enter  into  the  present 
investigation. 1 4 

14  C.  L.  Oberg,  "Acoustic  Instability  in  Propellant  Combustion,"  unpub¬ 
lished  Ph.D.  Thesis,  Department  of  Chemical  Engineering,  Univer¬ 
sity  of  Utah,  Salt  Lake  City,  Utah,  June  1965,  p.  72. 
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It  is  interesting  to  consider  the  attenuation  versus  pressure 
parametrically  described  by  the  ♦  ~<e  during  the  burn.  Such  a  graph 
corresponding  to  the  firing  represented  in  Fig.  7  is  given  in  Fig. 

8,  where  the  theoretical  equilibrium  attenuation  versus  pressure  is 
shown  as  a  dashed  line.  It  is  obvious  that  the  equilibrium  predic¬ 
tions  of  attenuation  do  not  follow  actual  results  during  initial  pres¬ 
surization  and  pressure  termination;  nevertheless,  there  exists  a 
period  of  time  during  which  the  measured  attenuation  agrees  with  the 
theory.  It  will  be  assumed  that  during  this  interval,  the  combustion 
gas  is  in  an  equilibrium  state  and  that  the  plasma  electrons  are  in 
thermal  equilibrium  with  the  gas.  It  then  appears  proper  to  ascribe 
the  non-correlation  points  of  Fig.  8  to  a  dynamic  state  of  plasma  for 
which  the  theory  does  not  apply.  The  attenuation  seems  to  lead  the 
pressure  rise  and  lag  the  pressure  decrease  at  the  end  of  the  bum, 
and  while  the  measured  attenuations  present  interesting  questions  per¬ 
taining  to  electron  temperature  and  gas  composition  during  ignition 
and  burnout,  the  developed  theory  cannot  be  used  to  predict  the  micro¬ 
wave  beam  interaction  for  these  conditions;  hence,  it  must  only  be 
applied  during  the  times  when  an  equilibrium  condition  may  be  safely 
assumed . 

The  graph  of  Fig.  8  shows  more  time  of  assumed  equilibrium  at 
the  end  of  the  burn  than  at  the  beginning.  It  seems  justified  to  ac¬ 
count  for  a  lot  of  the  relatively  long  time  between  ignition  and  equi¬ 
librium  in  the  thermal  mass  and  characteristics  of  the  heavy  motor 
casing. 
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Fig.  8.  Plotted  data  corresponding  to  Fig 


Figure  9  shows  experimental  attenuation  measurement  results 
for  a  5  percent  ?eeding  of  KCIO^.  The  representation  of  Fig.  9  and 
chat  of  the  subsequent  three  figures  follows  the  format  of  Fig.  8  in 
showing  the  time  evolution  of  attenuation  during  the  burn.  An  abbre- 
.iated  form  of  the  firing  shown  in  Pig.  8  is  seen  in  Fig.  9  on  the 
far  left.  Figure  10  depicts  experimental  data  for  the  5  percent  seed¬ 
ing  versus  equilibrium  temperature.  Figures  11  and  12  show  experimen¬ 
tal  data  for  a  2  percent  seeding  of  KC1Q.  and  correspond  to  Figs.  9 
and  10.  The  four  graphs  clearly  show  a  very  good  agreement  with  the 
theoretical  attenuations. 

it  will  be  noted  that  the  presented  experimental  results  agree 
with  the  theoretical  predictions  based  on  the  assumption  that  there 
exists  both  sources  and  sinks  of  electrons  in  the  plasma.  The  assump¬ 
tion  that  no  electron  sinks  exist  leads  to  predicted  attenuations  that 
are  an  rder  of  magnitude  higher  than  the  measured  results.  It  is 
then  obvious  that  the  chlorine,  acting  as  an  electron  sink,  plays  an 
important  role  in  determining  the  amount  of  attenuation  in  the  cham¬ 
ber. 

Th£  logical  conclusion  t  j  this  section  is  that  for  the  stated 
conditions,  the  theory  of  Chapter  II  does  represent  an  acceptable  mathe¬ 
matical  model  of  an  electromagnetic  plane,  wave  penetrating  a  homogen¬ 
eous  plasma.  The  practical  consequences  of  this  conclusion  may  now  be 
presented  to  indicate  some  of  the  uses  of  the  theory. 
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2  percent  potassium  perchlorate 


3.2  Consequences  of  the  Ver¬ 
ification  of  the  Theory 


Using  the  assumption  that  the  developed  model  is  correct,  many 
relationships  are  seen  to  exi?t  among  the  four  quantities;  attenuation, 
pressure,  temperature,  and  potassium  perchlorate  concentration.  In 
addition,  some  transient  conditions  may  be  analyzed  under  restricted 
circumstances  even  though  the  periods  of  ignition  stabilization  and 
burnout  represent  a  presently  undefined  stat°  of  plasma  and  cannot  be 
considered.  The  steady-state  relationships  are  considered  first. 

A.  Steady-State  Analysis 

When  the  four  quantities  noted  above  are  taken  as  variables,  the 
results  show  that  for  a  plasma  meeting  the  specified  requirements,  the 
determination  of  any  two  of  the  variables  will  establish  the  value  of 
the  other  two.  A  brief  description  along  with  a  possible  application 
will  be  given  of  each  situation. 

a.  Given  a  Measured  Attenuation  and  KC10/f  Concen¬ 
tration,  Find  the  Pressure  and  Temperature 

This  method  affords  a  verv  simple  way  to  determine  the  thermody¬ 
namic  state  of  the  gas  from  the  attenuation  alone  provided  the  KCIO^  is 
closely  controlled.  The  wave  generating  and  measuring  apparatus  is  the 
only  equipment  required  in  this  instance  allowing  this  system  tc  be  ad¬ 
vantageous  for  very  high  pressures  where  pressure  is  not  easily  meas¬ 
ured  directly.  This  method  would  be  useful  if: 

1.  A  sufficient  quantity  of  electrons  could  be  produced  at  the 
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2.  The  transmission  windows  could  be  made  to  withstand  the 
high  pressure. 


b.  Given  a  Measured  Attenuation  and  Pressure, 

Find  the  Temperature  and  Amount  of  KCIO^ 

These  conditions  constituted  the  basis  of  the  experimental  test 
unit  except  that  the  percent  KCIO^  was  known.  Since  the  attenuation, 
however,  was  found  to  be  raonctonic  with  respect  to  both  pressure  and 
KCIO^  concentration,  one  attenuation  measurement  for  a  given  pressure 
will  determine  both  the  temperature  and  KCIO^  in  the  propellant.  This 
choice  of  variables  would  be  useful  in  the  correlation  and  calibration 
of  other  temperature  measuring  devices  and  would  be  helpful  in  the 
study  of  the  thermal  equilibrium  characteristics  of  different  combus¬ 
tion  chamber  configurations.  A  plot  similar  to  Fig.  8  would  determine 
the  length  of  time  required  to  establish  an  equilibrium. 

The  determination  of  the  KC10.  concentration  could  be  used  as 
a  quality  control  check  in  the  production  processing  of  the  propellant. 
The  attenuation  can  correctly  indicate  the  Included  amount  of  potas¬ 
sium  or  any  other  ionizable  substance,  and  if  a  studied  constituent  of 
the  composition  could  be  made  proportional  to  a  reasonable  cmount  of 
electron  producing  material,  the  attenuation  measurements  from  the  bum 
of  a  ba'.ch  sample  would  indicate  the  amount  of  the  studied  component 
present .. 
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C.  GiVCS  «  MCuoul Gu  At tenUgu xun  auu  icaipetaiutCi 
Find  the  Pressure  and  KClO^  Concentration 

This  method  could  be  used  in  the  absence  of  a  pressure  trans¬ 
ducer  to  find  the  pressure  for  a  given  attenuation  if  the  temperature 
could  be  established  using  other  than  microwave  means.  The  amount  of 
KCIO^  is  also  established,  and  the  utility  of  this  information  is  as 
listed  in  b;  however,  the  difficulty  in  measuring  the  plasma  tempera¬ 
ture  would  seem  to  limit  this  type  of  system. 

d.  Given  a  Measured  Pressure  and  Temperature, 

Find  Attenuation  and  KCIO^  Concentration 

These  conditions  are  the  converse  of  tt.ose  in  a  and  predict  the 
behavior  of  the  wave  in  the  plasma  and  amount  of  ionized  material  for 
a  known  gas  state.  The  prediction  of  wave  behavior  could  be  of  impor¬ 
tance  in  designing  circuits  utilizing  the  plasma  as  an  element,  for 
when  the  power  loss  in  the  plasma  is  determined,  the  system  power  level 
can  be  optimized  with  respect  to  the  rest  of  the  system.  Unfortunately, 
the  assumption  of  a  given  temperature  suffers  the  disadvantages  listed 
in  c,  so  wave  behavior  is  best  determined  in  the  following  way  which 
assumes  no  known  temperature. 

e.  Given  a  Measured  Pressure  and  KCIO^  Concen¬ 
tration,  Find  Attenuation  and  Temperature 

For  a  controlled  amount  of  KCIO^  as  in  a,  only  the  pressure  is 
required  to  determine  wave  behavior.  The  temperature  is  also  estab¬ 
lished  in  this  method. 
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f.  Given  a  Measured  Temperature  and  KCIO^  Con¬ 
centration,  Find  Attenuation  and  Pressure 

For  systems  in  which  the  temperature  might  be  more  reliably 
known  than  pressure,  this  method  could  be  used  to  find  the  wave  at¬ 
tenuation  and  pressure  where  the  amount  of  KCIO^  is  known. 

As  a  summary  to  the  possibilities  presented,  it  will  be  seen 
that  there  are  three  ways  to  obtain  each  of  the  four  variables,  each 
way  containing  inherent  practical  considerations  yet  contributing  an 
avenue  of  solution.  Of  the  most  easily  seen  features  of  the  results 
is  that  the  attenuation  becomes  low  and  flat  for  high  pressures  and 
high  and  steep  for  low  pressures.  The  system  parameters  must  be  ad¬ 
justed  to  allow  sufficient  confidence  in  and  measurability  of  the  re¬ 
sulting  attenuation.  In  correspondence  to  this,  the  temperatures  at 
high  pressures  change  by  smaller  amounts  than  at  low  pressures  for  a 
given  pressure  shift,  and  the  resulting  attenuation  changes  could  pre¬ 
sent  a  problem  in  resolution. 

Another  consideration  implicit  in  the  calculation  is  thci  the 
attenuation  sharply  increases  as  the  electron  density  becomes  great 
enough  to  make  the  plasma  frequency  approach  the  incident  wave  fre¬ 
quency;  therefore,  while  four  original  variables  were  assumed  as  sys¬ 
tem  parameters,  the  wave  frequency  is  seen  to  require  optimization. 

B.  Transient  Plasma  Analysis 

To  ascertain  the  usefulness  of  the  theoretical  conclusions  in 
the  transient  plasma  condition,  the  sources  of  the  imposed  limitations 
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are  nrst  nxscusseu.  lniiidxxy }  iuclc  is  cite  cowpllcatcd  PcriCo  of 
events  subsequent  to  the  instant  of  ignition  of  the  surface  of  the 
propellant.  At  ignition,  the  surface  outgases  any  contaminants  that 
•cay  have  diffused  into  the  surface,  and  it  thus  represents  a  thin 
layer  of  unique  substance  because  of  this  diffusion  and  also  because 
of  chemical  reactions  of  the  surface  with  the  atmosphere. 

The  chamber  volume  at  ignition  is  filled  with  room  temperature 
air  which  atmosphere  then  becomes  mixed  with  increasing  quantities  of 
combustion  products  until  pressurization  initiates  nozzle  flow. 

The  heavy  chamber  structure  at  ignition  is  at  room  temperature 
and  represents  a  considerable  heat  sink  to  the  gas  molecules  which  soon 
rise  to  a  temperature  of  nearly  ten  times  room  temperature. 

The  complexity  of  the  three  initial  factors;  the  inhomogeneous 
propellant  surface,  the  rapidly  changing  chamber  gas  composition,  and 
the  large  thermal  mass  of  the  chamber;  easily  shows  why  the  attenuation 
measured  during  pressurization  cannot  be  described  by  the  theory.  It 
appears  that  the  order  of  stabilization  of  these  three  factors  in  time 
is  the  order  given  above,  respectively,  for  three  reasons: 

1.  •-  Ignition  of  the  highly  sensitive  surface  appears  to  be  almost 

instantaneous . 

2.  The  pressure  rise  time  from  10  percent  to  90  percent  of  peak 
value  has  exhibited  time  periods  of  the  order  of  one  second. 

3.  The  total  time  taken  for  the  measured  attenuation  to  approach 
the  theoretical  value  is  several  times  the  press  are  rise 
time. 
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It  can  then  be  assumed  that  after  ignition,  the  propellant  first 
becomes  a  homogeneous  burning  surface,  then  the  chamber  volume  becomes 
a  homogeneous  though  time-changing  atmosphere,  and  finally  the  thermal 
mass  of  the  case  becomes  a  constant  heat  transfer  medium.  During  the 
initial  stabilization  of  the  plasma  in  the  chamber,  the  attenuation 
results  from  a  series  of  nearly  totally  undefined  conditions  of  com¬ 
bustion,  atmosphere,  and  temperature,  and  might  be  expected  to  disa¬ 
gree  with  the  steady  case.  Indeed,  this  is  true,  as  Fig.  8  shows  the 
initial  attenuation  characterized  by  a  high  early  peak  and  a  greater 
signal  irregularity  than  during  any  other  time  of  the  burn. 

Conditions  at  the  time  of  propellant  burnout  are  of  a  similar 
complexity,  and  ambiguities  arise  from  the  transition  of  a  hot,  dense 
gas  to  a  stagnant  pool  of  spent  reactants  and  chamber  residues.  As 
combustion  terminates,  the  components  change  with  adherence  to  consid¬ 
erations  of  inst  taneous  nozzle-determined  pressure,  the  rapidly  cool¬ 
ing  chamber  walls,  and  the  reaction  rates  of  the  components.  These  con¬ 
ditions  again  make  it  impossible  to  predict  the  attenuation,  and  they 
explain  why  the  attenuation  quickly  diverges  from  the  theoretical  re¬ 
sult  once -combustion  has  ceased. 

The  conclusion  of  the  transient  analysis  limitations  show  that 
the  developed  model  cannot  be  used  in  the  analysis  of  large  changes  of 
state  in  the  chamber;  however,  those  fluxuations  occurring  when  an  equi¬ 
librium  condition  is  apparent, that  neither  alter  the  gas  composition 
much  nor  change  the  temperature  appreciably  may  be  analyzed.  One  such 
transient  event  allowing  analysis  is  pressure  variations  and  irregu- 
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larities.  During  several  experimental  burns,  fragments  momentarily 
clogging  in  the  nozzle  caused  sudden,  brief,  and  measurable  pressure 
rises  and  corresponding  instantaneous  dips  in  attenuation.  Unfortu¬ 
nately,  the  apparatus  was  not  designed  with  the  capability  of  chang¬ 
ing  the  chamber  pressure  during  burn,  so  this  feature  could  nor  be 
investigated;  nevertheless,  there  is  strong  indication  that  this 
method  could  be  used  to  analyze  nozzle  characteristics  and  degrada¬ 
tion  in  the  absence  of  a  pressure  transducer. 

Another  transient  analysis  that  could  be  performed  on  the 
plasma  is  the  time  analysis  of  different  regions  throughout  the  cham¬ 
ber.  The  periods  of  equilibrium  could  be  determined  as  in  Fig.  8  for 
various  cross  sections  in  the  plasma  provided  the  beam  could  traverse 
the  path  without  appreciable  diffraction.  Such  an  analysis  might 
facilitate  the  design  of  combustion  chambers  by  yielding  information 
about  the  thermal  equilibrium  time  lag  throughout  the  structure. 

This  brief  summary  of  possible  applications  of  the  theory  indi¬ 
cates  that  the  microwave  diagnostic  method  provides  powerful  tools  in 
the  evaluation  of  many  properties  of  the  entire  system  of  plasma  and 
container.  Once  the  parameters  of  beam  wavelength  and  ionizing  mate¬ 
rial  have  been  optimized  with  respect  to  the  resolution  of  the  measur¬ 
ing  equipment,  a  very  critical  analysis  of  the  rocket  motor  system  is 
afforded. 
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3. 3.  Conclusions  and  Contribution 
of  the  Investigation 


The  fundamental  conclusion  of  this  investigation  is  that  there 
exists  a  useable  relation  between  the  equilibrium  thermodynamic  state 
of  a  given  homogeneous,  combustion  chamber  plasma  and  the  attenuation 
of  an  incident  plane  electromagnetic  wave  of  wavelength  small  with 
respect  to  the  plasma  dimensions. 

The  consequences  of  this  conclusion  have  indicated  many  ways  in 
which  a  microwave  signal  can  be  used  to  analyze  and  optimize  the  opera-* 
cion  of  a  rocket  motor.  The  plotted  results  show  that  a  change  of  sev¬ 
eral  degrees  of  temperature  out  of  several  thousand  can  be  easily  de¬ 
tected  by  the  incident  beam;  moreover,  the  only  perturbation  of  the 
measurement  is  the  requirement  of  transmission  windows  in  the  chamber 
walls  and  a  measurable  number  of  free  electrons  in  the  combustion  gas. 
Some  of  the  advantages  of  the  microwave  analysis  of  a  chamber  are  the 
instantaneous  response  of  the  system,  the  ability  of  analysis  without 
other  transducers,  and  spacial  resolution  in  studying  localized  dis¬ 
turbances  in  the  chamber. 

The  basic  assumptions  of  the  method  require: 

1.  That  a  plane  wave  propagates  through  a  plasma  of  essentially 
negligible  reflection  coefficient. 

2.  That  the  electron  collision  phenomena  be  defined. 

3.  That  the  number  of  free  electrons  be  measurable. 

With  respect  to  the  first  requirement,  a  plane  wave  can  easily  be  ob¬ 
tained  it  proper  microwave  techniques  are  followed;  the  reflection  co~ 
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efficient  is  no  problen  for  a  gas  whose  plana  frequency  is  lower  than 
the  incident  wave  frequency.15  That  the  second  requirement  can  be  sat** 
lsfled  in  general  is  evidenced  by  the (ability  of  the  complex  collision 
frequency  analysis  (Eq.  1.12>  to  produce  verified  results  over  a  range 
of  pressure.  To  establish  a  nuaber  of  free  electrons  in  the  plasma  was 
the  reason  KCIO^  was  added  to  the  propellant.  This  addition  effectively 
swamped  the  random  generation  of  electrons,  but  the  large  amount  of 
chlorine  had  to  be  considered  as  a  source  of  considerable  electron  re¬ 
combinations.  Again,  the  validity  of  the  method  used  to  determine  the 
number  density  of  electrons  is  shown  by  the  verified  results  over  a 
wide  range  of  KCIO^  amounts.  It  can  therefore  be  concluded  that  the 
assumptions  of  this  technique  do  not  in  general  pose  a  problem  where 
combustion  plasmas  are  concerned. 

The  contribution  of  this  Investigation  is  the  manifestation  of 
the  usefulness  of  the  Incorporation  of  a  microwave  transmission  system 
into  the  evaluation  of  the  thermodynamic  properties  of  the  combustion 
plasma  of  s  solid  propellant  rocket  motor.  The  evidence  indicates  that 
the  microwave  method  has  a  definite  place  among  the  tools  of  rocket 
analysis. 


15  G.  Gal  and  W.  Gibson,  "Electromagnetic  Wave  Interaction  Tilth  Cylin¬ 
drical  Plasma,"  IEEE  Transactions  on  Antennas  and  Propagation,  Vol. 
AP-16,  July  1968,  pp.  468-475. 
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THE  ATTENUATION  PROGRAM  LISTING  AND 
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31  rCm  MAIN 

C  ATTENUATION  PROGRAM 

DlNLi&iw'.  HHC«iii»k>w^.<n>»onUil>(XHU)«X(U)»Mf>ai-)  >ATIN(11)» 

lATTNUBUlt  »TC11> 

REAL  M0FH2*  MOFH20,  uO»  N2»  MOFC02»MOFHCL»MOFCO.PMFK2*MMFK20* 

1  MMFN2f MMFC02*MMFHCl  ,K'!FCO,HOTuT,MMFK 

complex  p»c»r* ycurf < i; 

C  HEAD  IN  all  CONSTANTS 

UATA  uLANK»ZMaSS.BCLT2»ZCHAHG.PI»  AirOS.  FREW.  QVSH20. 
lQvSHCL.  UVC02»OCO»0»I?.OOVCO»OOVH2.  OCVN2/6.62517E-34. 
29.10B3E-31. 1.38044E-23.  l.to02tt“’l9»3.l4l596»6»0224E23»22«72:>E9» 
3b.92E**»l.335£-a»7.t6E-14»2#65E“20t8.92E-2u»2.08E-25»1.464!£-25» 
42.7t-25/ 

W=2.*PI*FREQ 

NSQ=tl**2 

CLOT  s2 • 99793E8 

CK=W/C‘.GT 

C=8./(3.*S0RT(PI)I 

Z»1A=6.0224E23 

EV=4.31 

07=0.05 

JDJS1 

C  REAP  in  THE  PERCENT  OF  POTASSIUM 

48  READ  50.PERCNT 
IFIPERCIITl  7»7*47 

47  00  25  i= 1.6 

C  RE AO  IN  PRESSURE  ANn  MCLE  FRACTION  OATa 

49  READ  50.  PR£SP.M0T0T*MMFH2.NmFH20.MMFN2»MMFc02»MMFHCL.MHFC0»MMFK 
HEAD  50. CL 

IFIPRESP:  7.7.51 
51  K=1 

PKtSsPRESP/14.7 

MOFHir  HMFHVMOTOT 

M0FH2C=HHFH2O/M0T0T 

MOFN2SMMFN2/MOTOT 

MOFC02=MMFC02/MOTOT 

MOFHCLSMMFHCL/MOTOT 

MOFCOsMMFCO/MOTOT 

Zlt=MMFK/HOTOT«ZHA 

ZUCL=CL/MOTOT«ZKA 

TEMPS2.E3 

C2=M0FH20*QVStl20«M0FHCL*0VSHCL 
C l=M0FCO2*OVCO2 
CZ=MOFCO*OCO+MOFH2*OH2 

CM=MOFH2*OuVH2+MOFN?*OOVN2*MOFCO*QOVCO 
1  IFIK-12)  69.2.2 

C  CALCULATE  THE  INTE6KAL.V  0AR=B-J0 
69  BETa=ZMASS/|2.*H0LT7*TEMPI 

RM0(K)SPHCS*AV06/(8?.05*TEMP»«1.E6 

TW=W*UETA 

T2=C2*RHC(K)»(S0K7 (OFTA* )**3 
71=C1«HH0<K)*UE?A 
TZ=CZ»RHO ( K ) *SORT I UFTA ) 

TM=CM*KmO(K) 

01=0. 

ni=o, 

HSU. 

u=o. 
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J=0 

50  Y=Y*UY 
JsJ+JOJ 

IP ( J-A5)  70.1ti0,10u 
70  0EN0MlS<?#*y)**2 
8L=B 
0L=0 

YP0LY=T2*TI*Y*TZ*Y*«;>*TM*Y*#s 
O£N0K2=YPOLr**2 
UCNOM=OL'NOH  i+OEl  iOM2 

OrC*Y«JOLY*Y*#b*EXP  I  -Y*#2 )  /DENOM*HETA 
0=C*St*0eTA*Y**6*rxP(-Y**2>/0EN0M*nETA 

CfiSO  .  5*  { JH.+8  >  *OY 
CDUr .  5*  C  OL+D i *0 Y 
(tl-UIlOU 
DIsDIiDUO 
GO  TO  30 
loo  ro(K)=tji 
DS(K)=01 

:  I0NIZ*T10N  AWO  PLASMA  FREQUENCY 

( “*CMAR6*EV/ <  HOLTZ*  T£MP ) ) 
f  rf  L^£XP  *  “ZCHAR°*3 .76/1 OOL  TZ*  TE MP  > » 

ZTEMPrT£MP*sai<T  { TEMP , 

f  52li«  2*  *PI  *Q* 1  u83* 1  •  >  z<>  •  ft25l  7)  **3*  1 .  £21 

ALPHAl-ZCONST*ZTEMrVZtXPK*i*F«2G 
ALPHa2=4.*2c0NST*ZTPMH»ZEXPCI  iI.E-20 

Al=ZK*l.£-20 

A2=ZNCL*l.E-2u 

Ps  CMPLX(ALPHAl+ALPHA2<A2,0.(i, 

Q=  CHPLX-.ALPHAl*IA2-Al+ALPilA2)#0.OI 
US  CMPLX(-ALPHA1»ALPHA2*A1,U.D) 

CALL  CUi?IC(P,Q,H,YCUHE; 

UO  500  LLS1.3 

LZ=LL 

I'*(REAI.IYCUBE(LZH)  hGu»bUU»Si;2 
500  CONTINUE 

PRINT  99 • TEMP 

"  x?Kj!-o,bX#,fW  S0LUT f 0NS  AKfc  POUND  FOR  A  TEMPERATURE  OF 

RMOtoory, 

WPtK)=0. 

ATTN'K)=0. 

ATTH0B(K>s(|# 

GO  TO  37 

502  RMOE(M=REAL(YCUBEU  7)  »*l.fc2rt*2CMARG 
XCM  IS  *  ME  NUMBER  DFNSITY  OF  ELECTRONS 
X<K|=REAL(YCUUE(LZ) >*1.E2u 
2I€.=RH0EIK) 

=2»JE*I.  602/9.  1 0B3*3b.  *Pl*  I  .£21 
WP(K)  la  THE  PLASMA  FREQUENCY 
WP  (  K  }  =SQRT  ( AtiS  <  WPSu  1 » 

T(K>=TEMP 

c‘iCUL,TE  T,c  *,reNi"™" 
«im  kT"  '  1  “ '  ’  *  ’ 2  ’ 

IF<TERM1+TERM2)  bCbrSlu.blO 


*»F6.1//} 
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5u5  Pm  NT  oOb.TEMP 

606  FL»RM*T!bX;*  ATTENUATION  72™  fS  NEGATIVE  FOR  a 


TEMPERATURE  OF 


lFti.l//) 

ATTN<K)=(j. 

ATTNOU(K)=0. 

GO  TO  37 

blO  ATTU{N)=SORTtTERMl*TFRM2>*CK/SOHT<2.) 

ATTS<OB  IS  UH  ATTENUATION  PfcH  2.20  INCH  ATTENUATION  PATH 
ATTUUb  C  N  )  =20  .  *  .  i3‘»2‘J4a«*2 . 54F~2*2 . 20  *  ATTN  CK  ) 

37  T<M=TEMP 

TtMP=TEHP+ioo. 

K=K*  1 
GO  TO  1 

2  PRINT  3,PRESP.P£KCNT 
PRINT  4 

PRINT  5 

PRINT  O,  tTlj) »RHO( J1 ,KHO£( J) .HBt J> .DD( J) ,X(  J) ,WP< J) ,ATTN{ J) .ATTNOR 

l(sl).J=l»ll  J 

25  PRINT  10 
GO  TO  4d 
7  CONTIIAJL 
10  FORMAT (1H1) 

50  FORMAT (FO.l.dFS.b) 

3  FORMAT C17X»»CALCULATFU  ATThNHAT ION  PARAMETERS  FOR  ’»F5.1»*  PSIA  AN 
1U  *»F4.1,»  PERCENT  POTASSIUM  PERCHLORATE •/✓/ I 

4  FORMAT <  IX,  • Ttftf** • •  GAS  *,*  CHARGE  ».*  V  BA 

ir  integration  *,*  electron* »•  plasma  •»» 

2ATTLNUATI0N  •> 

5  FORMAT (5X»  *  DENSITY  DENSITY  » ,9X, *B» * 13X. *D* »6X» • 

1  DtNSiTy  *.*  FREQUENCY  •**  NEPERS/MCTER  DB*  *//) 

6  FORMAT { lX#F7.1*E12.7»faElb.7#F10.3//> 

END 


Ul  FOR  CUBIC 

SUBROUTINE  CUBIC IP* O.H*Y) 

C  *ROUTl!lt  TO  SOLVE  A  CUBIC  EQUTION. 
COMPLEX  P,Q,R,Y(3>,A.B,V,P0,71,T 
UAl A  M/(-0.5.0.HB6u?S404J/ 

P0=P/3.,0 
A=Q/3.0-P0**2 
B=P0**3*  0 . b*  ?  K-P0*0 1 
ZI--3 

IF { ABS (RLAL  5 A) ) *ABS ( AiMAb ( A) » )  1.3,1 

1  T=CSQRT  *  A**3*U**2 ) 

Z1=-B*T 

IF(RLAL(BH  3,3.2 

2  Zl=-B-T 

3  IF;AIMAG<Z1> )5.4,5 

4  Z1=CMPLX(CBRT(REAL(71 ) ) ,U.U» 

GO  TO  u 

5  Zl=tCt)RT(Zl) 

6  NO  7  K=i , 3 
Y(K)=Z1-A/Z1-P0 

7  Z1=W»Z1 
RETURN 
ENU 

W  XQT  MAIN 


CUBIC  l 
CUBIC  2 
CUBIC  3 
CUBIC  •> 
CUBIC  5 
CUBIC  6 
CUBIC  7 
CUBIC  6 
CUBIC  9 
CUBIC  10 
CUBIC  11 
CUBIC  12 
CUBIC  13 
CUBIC  14 
CUBIC  15 
CUBIC  16 
CUBIC  17 
CUBIC  18 
CUBIC  19 
CUBIC  20 
CUBIC  21 
CUBIC  22 


I.  **»»*>:  /  . . . ‘11  •  . .  1  ""  *  . . 


COwnUSTiOli  MIASMA  COMPOSITION  DATA 
1.0  PENCE! IT  KCLl/4  0Y  WEIGHT  * 

CHAMUEK  TOTAL  MCLES  MOLES  MOlES  MOLES  MOLES  MOLES  MOLES  MOLES 
PRESSIVE  MOLES  OF  OF  Of  OF  OF  OF  OF  OF 


PSIA 

GAS 

M2 

H?0 

H2 

C02 

HCL 

CO 

K 

CL 

120. 

4.227 

.oioua 

l.3u ins 

•  3S9li7 

.30591 

.67232 

.91438 

.00905 

.01709 

luo. 

4.229 

.61*160 

1.30077 

<  3.69U5 

.30395 

.67110 

.91435 

.00005 

.01951 

8u. 

4.232 

.61137 

1.2994U 

.35903 

*3u6u3 

.66951 

.91426 

.90006 

.01990 

t>u. 

4.23S 

.61231 

1.29754 

.35699 

.30616 

.66733 

.51415 

.00006 

. 02209 

4o. 

4.2** 

.61367 

1. ?9472 

.35695 

.30642 

.66396 

.91391 

.00008 

.02546 

2u. 

4.2S1 

.61ol9 

1.26934 

.3586b 

.3070? 

.65751 

.91333 

.0001 

.03195 

2.0  PERCENT  KCL04  OY  WEIGHT  * 

CHAMUEK  TOTAL  MOLES  HOLES 
PRESSURE  MOLES  OF  OF 


MOLES  MOLES  MOLES  MOLES 
OF  OF  OF  OF 


MOLES  MOLES 
OF  OF 


PSIA 

GAS 

H2 

H20 

112 

120. 

4.2U9 

.  OO  _  7  .i 

1.29^55 

.3:n»82 

too. 

4.211 

.00^29 

f .29429 

.35480 

SHu. 

“.215 

.60500 

1.29293 

.35477 

u6u. 

4.217 

.00593 

1.29Ju7 

. 35474 

u4u. 

4.222 

.6072° 

1 .2662- 

.354b9 

02U. 

4.232 

.00979 

1.26207 

„354bl 

CO? 

HCL 

CO 

CL 

3**7 1 1 

.66403 

.91318 

.00310 

.01697 

3o7lb 

. S6282 

.91314 

.00011 

.01818 

3u724 

'66125 

.91307 

.00012 

.01976 

30737 

•b5909 

.91295 

.00013 

.02193 

30762 

.65577 

.91271 

.00015 

. 02527 

3u822 

.64937 

.91214 

.00020 

.03172 

COMBUSTION  PLASMA  COMPOSITION  OATA  CONT. 


•* 


5.0  PERCENT  KCL04  BY  WEIGHT  ♦ 


CHAMBER 

PRESSURE 

PSIA 

total 

MOLES 

GAS 

MOLES 

OF 

H2 

* - 

moles 

OF 

H20 

MOLES 

OF 

H2 

120. 

4*154 

.58472 

1.27582 

.34205 

100. 

4.15s 

.58527 

1.2747b 

.34203 

uBu. 

4.158 

•5859b 

1.27341 

.34200 

u6o. 

4.161 

.58687 

1 .27l5b 

.34197 

u4u. 

4*  Iu7 

.58820 

1.26875 

.34192 

u2u. 

4.177 

.59u65 

1.26337 

.34184 

MOLES 

OF 

C02 

MOLES 

OF 

HCL 

MOLES 

OF 

CO 

MOLES  MOLES 
OF  OF 

K  CL 

31033 

.63917 

.90946 

.f  325  .0166 

31087 

.63799 

.90943 

.00027  .01779 

31094 

•63648 

.90936 

.00030  .01933 

31107 

.63438 

.90925 

.00033  .02146 

31131 

.63117 

.90902 

.00039  .02473 

31189 

.62498 

.90846 

.00050  .03102 

10.0  PERCEPT  KCL04  BY  WEIGHT  * 

chamber 

PRESSURE 

PSIA 

total 

MOLES 

GAS 

MOLES 

OF 

H2 

MOLES 

OF 

H20 

MOLES 

UF 

H2 

MOLES 

OF 

C02 

moles 

OF 

HCL 

MOLES 

OF 

CO 

MOLES  MOLES 
OF  OF 

K  CL 

120. 

4.002 

.55333 

1,24294 

.32077 

.31732 

.59779 

.90297 

.00952  .01598 

100. 

4.004 

.55387 

1.24109 

.32075 

.31736 

.59668 

.90295 

.00056  .01712 

0Bu. 

4.006 

.55453 

l .24054 

.32073 

.’1743 

.59525 

.90269 

.00061  .01861 

U6u. 

4.  Qo9 

.55540 

1.2387g 

.32069 

. 31754 

.59327 

.90278 

.00069  .02065 

040. 

4.075 

.55oo8 

1 .2359u 

.32064 

.31777 

.59023 

.90256 

.09080  .02379 

U2U. 

4.055 

.559o6 

1.23054 

.32056 

.31831 

.58440 

,902p4 

.00104  .02983 
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I.  AM’RAC  T 


The  analysis  of  a  plane  wave  propagating  through  a  homogeneous  plasma  region  of 
dimensions  large  compared  to  the  signal  wavelength  has  shown  that  the  wave  attenuatioi 
can  be  related  to  physical  parameters  of  the  plasma  such  as  electron  density, 
temperature,  and  pressure.  In  this  report,  the  plasma  in  the  combustion  chamber  of  a 
small  solid  propellant  rocket  motor  is  used  to  investigate  the  feasbility  of 
determining  the  temperature  of  the  plasma  by  using  microwave  attenuation  data. 
Potassium  perchlorate  is  substituted  for  some  of  the  ammonium  perchlorate  ir  the 
propellant  to  give  an  electron  density  that  produces  an  easily  measurable  ictenuation 
Jhe  theoretical  results  show  that  for  a  known  propellant  composition  and  potassium 
perchlorate  concentration,  both  the  plasma  temperature  and  chamber  pressure  at 
equilibrium  may  be  found  from  one  attenuation  measurement.  Also,  for  a  known  thermo¬ 
dynamic  state  of  the  plasma,  the  electron  density  may  be  found  from  the  attenuation 
data.  The  method  of  microwave  analysis  Si  shown  to  be  verv  useful  in  the  study  of 
transient  phenomena  in  that  no  transducers  are  required;  hence,  the  time  response 
of  the  measurement  is  limited  only  by  the  microwave  detection  circuitry.  This 
technique  may  be  used  to  study  localized  distrubances  within  the  plasma  region. 
Experim#  ital  data  are  presented  to  support  the  theory  for  an  aranonium  perchlorate-  * 
polyfcu  dience  acrylic  acivi  propellant  which  was  burned  at  chamber  pressures  between 
20  and  80  psia  corresponding  to  a  temperature  range  of  2550°K  to  2600°K.  The  data 
show  that  temperature  changes  of  less  than  ten  degrees  can  easily  be  detected  by 
microwave  attenuation  measurements  when  the  propellant  is  seeded  with  2-10  percent 
potassium  perchlorate.  Seeding  otherpropellants  with  alkali  metal  compounds  is  expect 
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to  yield  comparable  results  ever  the  combustion 
range  in  temperature  and  pressure  of  the  other  pre 
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